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Abstract. It is well known that magnetic fields constrain motions of charged particles, impeding 
the diffusion of charged particles perpendicular to magnetic field direction. This modification of 
transport processes is of vital importance for a wide variety of astrophysical processes including 
cosmic ray transport, transfer of heavy elements in the interstellar medium, star formation etc. 
Dealing with these processes one should keep in mind that, in realistic astrophysical conditions, 
magnetized fluids are turbulent. In this review we single out a particular transport process, namely, 
heat transfer and consider how it occurs in the presence of the magnetized turbulence. We show that 
the ability of magnetic field lines to constantly change topology and connectivity is at the heart of the 
correct description of the 3D magnetic field stochasticity in turbulent fluids. This ability is ensured 
by fast magnetic reconnection in turbulent fluids and puts forward the concept of reconnection 
diffusion at the core of the physical picture of heat transfer in astrophysical plasmas. Appealing to 
reconnection diffusion we describe the ability of plasma to diffuse between different magnetized 
eddies explaining the advection of the heat by turbulence. Adopting the structure of magnetic 
field that follows from the modern understanding of MHD turbulence, we also discuss thermal 
conductivity that arises as electrons stream along stochastic magnetic field lines. We compare the 
effective heat transport that arise from the two processes and conclude that, in many astrophysically- 
motivated cases, eddy advection of heat dominates. Finally, we discuss the concepts of sub and 
superdiffusion and show that the subdiffusion requires rather restrictive settings. At the same time, 
accelerated diffusion or superdiffusion of heat perpendicular to the mean magnetic field direction is 
possible on the scales less than the injection scale of the turbulence. 
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MAIN IDEA AND STRUCTURE OF THE REVIEW 

Heat transfer in turbulent magnetized plasma is an important astrophysical problem 
which is relevant to the wide variety of circumstancies from mixing layers in the Local 
Bubble (see Smith & Cox 2001) and Milky way (Begelman & Fabian 1990) to cooling 
flows in intracluster medium (ICM) (Fabian 1994). The latter problem has been sub- 
jected to particular scrutiny as observations do not support the evidence for the cool gas 
(see Fabian et al. 2001). This is suggestive of the existence of heating that replenishes 
the energy lost via X-ray emission. Heat transfer from hot outer regions is an important 
process to consider in this context. 

It is well known that magnetic fields can suppress thermal conduction perpendicular 
to their direction. However, this is true for laminar magnetic field, while astrophysical 
plasmas are genetically turbulent (see Armstrong et al 1994, Chepurnov & Lazarian 
2010). The issue of heat transfer in realistic turbulent magnetic fields has been long 
debated (see Bakunin 2005 and references therein). 



Below we argue that turbulence changes the very nature of the process of heat trans- 
fer. To understand the differences between laminar and turbulent cases one should con- 
sider both motion of charged particles along turbulent magnetic fields and turbulent 
motions of magnetized plasma that also transfer heat. The description of both processes 
require the knowledge of the dynamics of magnetic field lines and the structure of the 
magnetic field lines in turbulent flows. The answers to these questions are provided by 
the theories of magnetic reconnection and magnetic turbulence. To provide the quanti- 
tative estimates of the heat transfer the review addresses both theories, discussing the 
generic process of reconnection diffusion which describes the diffusion induced by the 
action of turbulent motions in the presence of reconnection. We stress the fundamental 
nature of the process which apart from heat transfer is also important e.g. for removing 
magnetic field in star formation process (Lazarian 2005). 

In §2 we discuss the omnipresence of turbulence in astrophysical fluids, introduce 
major ideas of MHD turbulence theory and turbulent magnetic reconnection in §3 
and §4, respectively, relate the concept of reconnection diffusion to the processes of 
heat transfer in magnetized plasmas in §5. We provide detailed discussion of heat 
conductivity via streaming electrons in §6, consider heat advection by turbulent eddies 
in §7, and compare the efficiencies of the latter two processes in §8. Finally, we discuss 
heat transfer on scales smaller than the turbulence injection scale in §9 and provide final 
remarks in §10. 

MAGNETIZED TURBULENT ASTROPHYSICAL MEDIA 

Astrophysical plasmas are known to be magnetized and turbulent. Magnetization of 
these fluids most frequently arises from the dynamo action to which turbulence is an 
essential component (see Schekochihin et al. 2007). In fact, it has been shown that 
turbulence in weakly magnetized conducting fluid converts about ten percent of the 
energy of the cascade into the magnetic field (see Cho et al. 2009). This fraction does 
not depend on the original magnetization and therefore magnetic fields will come to 
equipartition with the turbulent motions in about 10 eddy turnover times. 

We deal with magnetohydrodynamic (MHD) turbulence which provides a correct 
fluid-type description of plasma turbulence at large scales 1 . Astrophysical turbulence is 
a direct consequence of large scale fluid motions experiencing low friction. This quantity 
is described by Reynolds number Re = LV/v, where L is the scale of fluid motions, V 
is the velocity at this scale and v is fluid viscosity. The Reynolds numbers are typically 
very large in astrophysical flows as the scales are large. As magnetic fields decrease 
the viscosity for the plasma motion perpendicular to their direction, Re numbers get 
really astronomically large. For instance, Re numbers of 10 10 are very common for 
astrophysical flow. For so large Re the inner degrees of fluid motion get excited and 
a complex pattern of motion develops. 



It is possible to show that in terms magnetic field wandering that is important, as we see below, for heat 
transfer the MHD description is valid in collisionless regime of magnetized plasmas (Eyink, Lazarian & 
Vishniac(2011). 



The drivers of turbulence, e.g. supernovae explosions in the interstellar medium, 
inject energy at large scales and then the energy cascades down to small scales through a 
hierarchy of eddies spanning up over the entire inertial range. The famous Kolmogorov 
picture (Kolmogorov 1941) corresponds to hydrodynamic turbulence, but, as we discuss 
further, a qualitatively similar turbulence also develops in magnetized fluids/plasmas. 

Simulations of interstellar medium, accretion disks and other astrophysical environ- 
ments also produce turbulent picture, provided that the simulations are not dominated 
by numerical viscosity. The latter requirement is, as we see below, is very important for 
the correct reproduction of the astrophysical reality with computers. 

The definitive confirmation of turbulence presence comes from observations, e.g. 
observations of electron density fluctuations in the interstellar medium, which produce 
a so-called Big Power Law in the Sky (Armstrong et al. 1994, Chepurnov & Lazarian 
2010), with the spectral index coinciding with the Kolmogorov one. A more direct piece 
of evidence comes from the observations of spectral lines. Apart from showing non- 
thermal Doppler broadening, they also reveal spectra of supersonic turbulent velocity 
fluctuations when analyzed with techniques like Velocity Channel Analysis (VCA) 
of Velocity Coordinate Spectrum (VCS) developed (see Lazarian & Pogosyan 2000, 
2004, 2006, 2008) and applied to the observational data (see Padoan et al. 2004, 2009, 
Chepurnov et al. 2010) rather recently. 

All in all, the discussion above was aimed at conveying the message that the 
turbulent state of magnetized astrophysical fluids is a rule and therefore the discussion 
of any properties of astrophysical systems should take this state into account. We shall 
show below that both magnetic reconnection and heat transfer in magnetized fluids are 
radically changed by turbulence. 

STRONG AND WEAK ALFVENIC TURBULENCE 

For the purposes of heat transfer, Alfvenic perturbations are most important. Numerical 
studies in Cho & Lazarian (2002, 2003) showed that the Alfvenic turbulence develops 
an independent cascade which is marginally affected by the fluid compressibility. This 
observation corresponds to theoretical expectations of the Goldreich & Sridhar (1995) 
theory that we briefly describe below (see also Lithwick & Goldreich 2001). In this 
respect we note that the MHD approximation is widely used to describe the actual 
magnetized plasma turbulence over scales that are much larger than both the mean 
free path of the particles and their Larmor radius (see Kulsrud 2004 and ref. therein). 
More generally, the most important incompressible Alfenic part of the plasma motions 
can described by MHD even below the mean free path (see Eyink et al. 2011 and ref. 
therein). 

While having a long history of ideas, the theory of MHD turbulence has become 
testable recently due to the advent numerical simulations (see Biskamp 2003) which con- 
firm (see Cho & Lazarian 2005 and ref. therein) the prediction of magnetized Alfvenic 
eddies being elongated in the direction of magnetic field (see Shebalin, Matthaeus & 
Montgomery 1983, Higdon 1984) and provided results consistent with the quantita- 
tive relations for the degree of eddy elongation obtained in Goldreich & Sridhar (1995, 
henceforth GS95). 



The hydrodynamic counterpart of the MHD turbulence theory is the famous Kol- 
mogorov theory of turbulence. In that theory, energy is injected at large scales, creating 
large eddies which correspond to large Re numbers and therefore do not dissipate en- 
ergy through viscosity 2 but transfer energy to smaller eddies. The process continues till 
the cascade reaches the eddies that are small enough to dissipate energy over an eddy 
turnover time. In the absence of compressibility the hydrodynamic cascade of energy is 
~ vj/T casc j = const, where v/ is the velocity at the scale / and the cascading time for the 

eddies of size / is x cas k.i ~ l/vi. From this the well known relation v/ ~ Z 1 / 3 follows. 

Modern MHD turbulence theory can also be understood in terms of eddies. However, 
in the presence of dynamically important magnetic field, eddies cannot be isotropic. 
Any motions bending magnetic field should induce a back-reaction and Alfven waves 
propagating along the magnetic field. At the same time, one can imagine eddies mixing 
magnetic field lines perpendicular to the direction of magnetic field. For the latter 
eddies the original Kolmogorov treatment is applicable resulting perpendicular motions 
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scaling as v// ± 7 , where lj_ denotes scales measured perpendicular to magnetic field and 
correspond to the perpendicular size of the eddy. These mixing motions induce Alfven 
waves which determine the parallel size of the magnetized eddy. The key stone of the 
GS95 theory is critical balance, i.e. the equality of the eddy turnover time l±/vi and 
the period of the corresponding Alfven wave ~ 1\\/Va, where /y is the parallel eddy 
scale and Va is the Alfven velocity. Making use of the earlier expression for v/ one can 
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easily obtain ~ l^ , which reflects the tendency of eddies to become more and more 
elongated as energy cascades to smaller scales. 

While the arguments above are far from being rigorous they correctly reproduce the 
basic scalings of magnetized turbulence when the velocity equal to Va at the injection 
scale L. The most serious argument against the picture is the ability of eddies to perform 
mixing motions perpendicular to magnetic field. We shall address this issue in §3 but 
for now we just mention in passing that strongly non-linear turbulence does not usually 
allow the exact derivations. It is numerical experiments that proved the above scalings 
for incompressible MHD turbulence (Cho & Vishniac 2000, Maron & Goldreich 2001, 
Cho, Lazarian & Vishniac 2002) and for the Alfvenic component of the compressible 
MHD turbulence (Cho & Lazarian 2002, 2003, Kowal & Lazarian 2010). 

It is important to stress that the scales lj_ and /y are measured in respect to the system 
of reference related to the direction of the local magnetic field "seen" by the eddy. This 
notion was not present in the original formulation of the GS95 theory and was added in 
Lazarian & Vishniac (1999) (see also Cho & Vishniac 2000, Maron & Goldreich 2001, 
Cho et al. 2002). In terms of mixing motions that we mentioned above it is rather obvious 
that the free Kolmogorov-type mixing is possible only in respect to the local magnetic 
field of the eddy rather than the mean magnetic field of the flow. 

GS95 theory assumes the isotropic injection of energy at scale L and the injection 
velocity equal to the Alfven velocity in the fluid Va, i.e. the Alfven Mach number 



2 Reynolds number Re = LV /v = (V /L)/(v/L 2 ) which is the ratio of the eddy turnover rate T edd = 

V/L and the viscous dissipation rate Tjj' = 77 /L 2 . Therefore large Re correspond to negligible viscous 
dissipation of large eddies over the cascading time T casc which is equal to l e ddy in Kolmogorov turbulence. 



Ma = (SV /Va) = 1. This model can be easily generalized for both Ma < 1 and Ma > 1 
at the injection (see Lazarian & Vishniac 1999 and Lazarian 2006, respectively). Indeed, 
if Ma > 1, instead of the driving scale L for one can use another scale, namely Ia, which 
is the scale at which the turbulent velocity gets equal to Va. For Ma 3> 1 magnetic fields 
are not dynamically important at the largest scales and the turbulence at those scales 
follows the isotropic Kolmogorov cascade v\ ~ Z 1 / 3 over the range of scales [L, Ia]. This 
provides Ia ~ LM A 3 . If Ma < 1, the turbulence obeys GS95 scaling (also called "strong" 
MHD turbulence) not from the scale L, but from a smaller scale l trans ~ LM^ (Lazarian 
& Vishniac 1999), while in the range [LJ trans ] the turbulence is "weak". 

The properties of weak and strong turbulence are rather different. The weak tur- 
bulence is wave-like turbulence with wave packets undergoing many collisions before 
transferring energy to small scales 3 . On the contrary, the strong turbulence is eddy-like 
with cascading happening similar to Kolmogorov turbulence within roughly an eddy 
turnover time. One also should remember that the notion "strong" should not be associ- 
ated with the amplitude of turbulent motions, but only with the strength of the non-linear 
interaction. As the weak turbulence evolves, the interactions of wave packets increases 
as the ratio of the parallel to perpendicular scales of the packets increases making the 
turbulence strong. In this case, the amplitude of the perturbations may be very small. 

While there ongoing debates whether the original GS95 theory should be modified to 
better describe MHD turbulence, we believe that, first of all, we do not have compelling 
evidence that GS95 is not adequate 4 . Moreover, the proposed additions to the GS95 
model do not change the nature of the physical processes that we present below. 

The quantitative picture of astrophysical turbulence sketched in this section gives 
us a way to proceed with the quantitative description of key processes necessary to 
describe heat transfer. The interaction of fundamental MHD modes within the cascade 
of compressible magnetized turbulence is described in Cho & Lazarian (2005), but this 
interaction is not so important for the processes of heat transfer that we discuss below. 

MAGNETIC RECONNECTION OF TURBULENT MAGNETIC 

FLUX 

Magnetic reconnection is a fundamental process that violates magnetic flux being 
frozen in within highly conductive fluids. Intuitively one may expect that magnetic fields 
in turbulent fluids cannot be perfectly frozen in. Theory that we describe below provide 
quantitative estimates of the violation of frozen in condition within turbulent fluids. 

We would like to stress that the we are discussing the case of dynamically important 
magnetic field, including the case of weakly turbulent magnetic field. The case of weak 
magnetic field which can be easily stretched and bended by turbulence at any scale up 



3 Weak turbulence, unlike the strong one, allows an exact analytical treatment (Gaultier et al. 2002). 

4 Recent work by Beresnyak & Lazarian (2010) shows that present day numerical simulations are unable 
to reveal the actual inertial range of MHD turbulence making the discussions of the discrepancies of 
the numerically measured spectrum and the GS95 predictions rather premature. In addition, new higher 
resolution simulations by Beresnyak (201 1) reveal the predicted —5/3 spectral slope. 
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FIGURE 1. Upper panel: Sweet-Parker reconnection. A is limited by resistivity and small. Middle 
panel: reconnection according to LV99 model. A is determined by turbulent field wandering and can be 
large. Lower panel: magnetic field reconnect over small scales. From Lazarian, Vishniac & Cho (2004). 

to the dissipation one is rather trivial and of little astrophysical significance 5 . At the 
same time, at sufficiently small scales magnetic fields get dynamically important even 
for superAlfvenic turbulence. 

Within the picture of eddies mixing perpendicular to the local magnetic field that 
we provided in the previous section, it is suggestive that magnetized eddies can provide 
turbulent advection of heat similar to the ordinary hydrodynamic eddies. This is rather 
counter-intuitive notion in view of the well-entrenched idea of flux being frozen in 
astrophysical fluids. As it is explained in Eyink et al. (201 1) the frozen-in condition is not 
a good approximation for the turbulent fluids 6 . The violation of the perfect frozenness 
of the magnetic field in plasmas also follows from LV99 model of reconnection (see 
discussion in Vishniac & Lazarian 1999). 

A picture of two flux tubes of different directions which get into contact in 3D space 



5 In the case of dynamically unimportant field, the magnetic dissipation and reconnection happens on the 
scales of the Ohmic diffusion scale and the effects of magnetic field on the turbulent cascade are negligible. 
However, turbulent motions transfer an appreciable portion of the cascading energy into magnetic energy 
(see Cho et al. 2010). As a result, the state of intensive turbulence with negligible magnetic field is short- 
lived. 

6 Formal mathematical arguments on how and why the frozen-in condition fails may be found in Eyink 
(2011). 



is the generic framework to describe magnetic reconnection. The upper panel of Figure 1 
illustrates why reconnection is so slow in the textbook Sweet-Parker model. Indeed, the 
model considers magnetic fields that are laminar and therefore the frozen-in condition 
for magnetic field is violated only over a thin layer dominated by plasma resistivity. The 
scales over which the resistive diffusion is important are microscopic and therefore the 
layer is very thin, i.e. A <C L x , where L x is the scale at which magnetic flux tubes come 
into contact. The latter is of the order of the diameter of the flux tubes and typically very 
large for astrophysical conditions. During the process of magnetic reconnection all the 
plasma and the shared magnetic flux 7 arriving over an astrophysical scale L x should be 
ejected through a microscopic slot of thickness A. As the ejection velocity of magnetized 
plasmas is limited by Alfven velocity Va, this automatically means that the velocity in 
the vertical direction, which is reconnection velocity, is much less than Va- 

The LV99 model generalizes the Sweet-Parker one by accounting for the existence 
of magnetic field line stochasticity (Figure 1 (lower panels)). The depicted turbulence 
is sub-Alfvenic with relatively small fluctuations of the magnetic field. At the same 
time turbulence induces magnetic field wandering. This wandering was quantified in 
LV99 and it depends on the intensity of turbulence. The vertical extend of wandering 
of magnetic field lines that at any point get into contact with the field of the other flux 
tube was identified in LV99 with the width of the outflow region. Note, that magnetic 
field wandering is a characteristic feature of magnetized turbulence in 3D. Therefore, 
generically in turbulent reconnection the outflow is no more constrained by the narrow 
resistive layer, but takes place through a much wider area A defined by wandering 
magnetic field lines. The extend of field wandering determines the reconnection velocity 
in LV99 model. 

An important consequence of the LV99 reconnection is that as turbulence amplitude 
increases, the outflow region and therefore reconnection rate also increases, which 
entails the ability of reconnection to change its rate depending on the level of turbulence. 
The latter is important both for understanding the dynamics of magnetic field in turbulent 
flow and for explaining flaring reconnection events, e.g. solar flares. 

We should note that the magnetic field wandering is mostly due to Alfvenic tur- 
bulence. To describe the field wondering for weakly turbulent case LV99 extended the 
GS95 model for a subAlfvenic case. The same field wandering 8 , as we discuss later, is 
important for heat transfer by electrons streaming along magnetic field lines. 

The predictions of the turbulent reconnection rates in LV99 were successfully tested 
3D numerical simualtions in Kowal et al. (2009) (see also Lazarian et al. 2010 for an 
example of higher resolution runs). This testing provided stimulated work on the theory 
applications, e.g. its implication for heat transfer. One should keep in mind that the LV 



Figure 1 presents only a cross section of the 3D reconnection layer. A shared component of magnetic 
field is going to be present in the generic 3D configurations of reconnecting magnetic flux tubes. 
8 As discussed in LV99 and in more details in Eyink et al. (201 1) the magnetic field wandering, turbulence 
and magnetic reconnection are very tightly related concepts. Without magnetic reconnection, properties 
of magnetic turbulence and magnetic field wandering would be very different. For instance, in the absence 
of fast reconnection, the formation of magnetic knots arising if magnetic fields were not able to reconnect 
would destroy the self-similar cascade of Alfvenic turbulence. The rates predicted by LV99 are exactly 
the rates required to make Goldreich-Sridhar model of turbulence self-consistent. 



model assumes that the magnetic field flux tubes can come at arbitrary angle, which 
corresponds to the existence of shared or guide field within the reconnection layer 9 . 

Alternative models of magnetic reconnection appeal to different physics to overcome 
the constraint of the Sweet-Parker model. In the Petcheck (1964) model of reconnection 
the reconnection layer opens up to enable the outflow which thickness does not depend 
on resistivity. To realize this idea inhomogeneous resistivity, e.g. anomalous resisitivity 
associated with plasma effects, is required (see Shay & Drake 1998). However, for 
turbulent plasmas, the effects arising from modifying the local reconnection events 
by introducing anomalous resistivity are negligible as confirmed e.g. in Kowal et al. 
(2009). Other effects, e.g. formation and ejection of plasmoids (see Shibata & Tanuma 
2001, Lorreiro et al. 2008) which may be important for initially laminar environments 
are not likely to play the dominant role in turbulent plasmas either. Therefore in what 
follows dealing with turbulent transfer of hear we shall appeal to the LV99 model of 
reconnection. 

RECONNECTION DIFFUSION AND HEAT TRANSFER 

In the absence of the frozen-in condition in turbulent fluids one can talk about recon- 
nection diffusion in magnetized turbulent astrophysical plasmas. The concept of recon- 
nection diffusion is based on LV99 model and was first discussed in Lazarian (2005) in 
terms of star formation 10 . However, reconnection diffusion is a much broader concept 
applicable to different astrophysical processes, including heat transfer in magnetized 
plasmas. In what follows we shall discuss several processes that enable heat transfer 
perpendicular to the mean magnetic field in the flow. 

The picture frequently presented in textbooks may be rather misleading. Indeed, 
it is widely assumed that magnetic field lines always preserve their identify in highly 
conductive plasmas even in turbulent flows. In this situation the diffusion of charged 
particles perpendicular to magnetic field lines is very restricted. For instance, the mass 
loading of magnetic field lines does not change to a high degree and density and 
magnetic field compressions follow each other. All these assumptions are violated in 
the presence of reconnection diffusion. 

We shall first illustrate the reconnection diffusion process showing how it allows 
plasma to move perpendicular to the mean inhomogeneous magnetic field (see Figure 
2). Magnetic flux tubes with entrained plasmas intersect each other at an angle and due to 
reconnection the identity of magnetic field lines change. Before the reconnection plasma 
pressure Ppi aS ma m the tubes is different, but the total pressure P p i aS ma + Pmagn is the 
same for two tubes. After reconnection takes place, plasma streams along newly formed 
magnetic field lines to equalize the pressure along two new flux tubes. The diffusion of 



y The model in LV99 is three dimensional and it is not clear to what extend it can be applied to 2D 
turbulence (see discussion in ELV1 1 and references therein). However, the cases of pure 2D reconnection 
and 2D turbulence are of little practical importance. 

10 Indeed, the issue of flux being conserved within the cloud presents a problem for collapse of clouds 
with strong magnetic field. These clouds also called subcritical were believed to evolve with the rates 
determined by the relative drift of neutrals and ions, i.e. the ambipolar diffusion rate. 



Before reconnection yz view After reconnection 




FIGURE 2. Diffusion of plasma in inhomogeneous magnetic field. 3D magnetic flux tubes get into 
contact and after reconnection plasma streams along magnetic field lines. Right panel: XY projection 
before reconnection, upper panel shows that the flux tubes are at angle in X-Z plane. Left Panel: after 
reconnection. 



plasmas and magnetic field takes place. The effect of this process is to make magnetic 
field and plasmas more homogeneously distributed in the absence of the external fields 1 1 . 
In terms of heat transfer, the process mixes up plasma at different temperatures if the 
temperatures of plasma volumes along different magnetic flux tubes were different. 

If turbulence had only one scale of motions its action illustrated by Figure 2 would 
create every flux tube columns of hot and cold gas exchanging heat with each other 
through the diffusion of charged particles along magnetic field lines. This is not the 
case, however, for a turbulence with an extended inertial cascade. Such a turbulence 
would induce mixing depicted in Figure 2 on every scale, mixing plasma at smaller and 
smaller scales. 

When plasma pressure along magnetic field flux tubes is the same, the connection of 
flux tubes which takes place in turbulent media as shown in Figure 3 is still important 
for heat transfer. The reconnected flux tubes illustrate the formation of the wandering 
magnetic field lines along which electron and ions can diffuse transporting heat. For the 



If this process acts in the presence of gravity, as this is the case of star formation, the heavy fluid 
(plasma) will tend to get to the gravitating center changing the mass to flux ratio, which is important to 
star formation processes. In other words, reconnection diffusion can do the job that is usually associated 
with the action of ambipolar diffusion (see numerical simulations in Santos de Lima et al. (2010). 
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FIGURE 3. Exchange of plasma between magnetic eddies. Eddies carrying magnetic flux tubes interact 
through reconnection of the magnetic field lines belonging to two different eddies. This enables the 
exchange of matter between eddies and induces a sort of turbulent diffusivity of matter and magnetic 
field. 



sake of simplicity, we shall assume that electrons and ions have the same temperature. In 
this situation, the transfer of heat by ions is negligible and for the rest of the presentation 
we shall talk about the transport of heat by electrons moving along wandering field 
lines 12 . 

Consider the above process of reconnection diffusion in more detail. The eddies 
1 and 2 interact through the reconnection of the magnetic flux tubes associated with 
eddies. LV99 model shows that in turbulent flows reconnection happens within one 
eddy turnover time, thus ensuring that magnetic field does not prevent free mixing 
motions of fluid perpendicular to the local direction of magnetic field. As a result of 
reconnection, the tube \ 1ow \\ u p transforms into 2 low 12 up and a tube 2 low 22 up transforms 
into \ 1 ™2V'P. If eddy 1 was associated with hotter plasmas and eddy 2 with colder 
plasmas, then the newly formed magnetic flux tubes will have both patches of hot and 



This is true provided that the current of diffusing hot electrons is compensated by the current of 
oppositely moving cold electrons, the diffusivity of electrons along wandering magnetic field lines is 
dominant compared with the diffusivity and heat transfer by protons and heavier ions. If there is no 
compensating current, electrons and ions are coupled by electric field and have to diffuse along wandering 
magnetic fields together and at the same rate. This could be the case of diffusion of plasmas into neutral 
gas. However, we do not discuss these complications here 
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FIGURE 4. Heat diffusion depends on the scale of the hot spot. Different regimes emerge depending on 
the relation of the hot spot to the sizes of maximal and minimal eddies present in the turbulence cascade. 
Mean magnetic field B is directed perpendicular to the plane of the drawing. Eddies perpendicular to 
magnetic field lines correspond to Alfvenic turbulence. The plots illustrate heat diffusion for different 
regimes. Upper plot corresponds to the heat spot being less than the minimal size of turbulent eddies; 
Middle plot corresponds to the heat spot being less than the damping scale of turbulence; Lower plot 
corresponds to the heat spot size within the inertial range of turbulent motions. 



cold plasmas. For the hierarchy of eddies the shedding of entrained plasmas into hot 
and cold patches along the same magnetic field lines allows electron conductivity to 
remove the gradients, conducting heat. This is the process of turbulent advection of heat 
in magnetized plasmas. 

The difference between the processes depicted in Figures 2 and 3 is due to the 
fact that the process in Figure 2 is limited by the thermal velocity of particles, while 
the process in Figure 3 depends upon the velocity of turbulent eddies only. In actual 
plasmas in the presence of temperature gradients plasmas along different elementary 
flux tubes will have different temperature and therefore two processes will take place 
simultaneously. 

Whether the motion of electrons along wandering magnetic field lines or the dy- 
namical mixing induced by turbulence is more important depends on the ratio of eddy 
velocity to the sonic one, the ratio of the turbulent motion scale to the mean free path 
of electrons and the degree of plasma magnetization. Strong magnetization both limits 
the efficiency of turbulent mixing perpendicular to magnetic field lines and the extent to 
which plasma streaming along magnetic field lines moves perpendicular to the direction 



of the mean field. However, but reduction of heat transfer efficiency is different for the 
two processes. We provide quantitative treatment of these processes in the next section. 

An interesting example of practical interest is related to the diffusion of heat from a 
hot spot. This case of reconnection diffusion is illustrated by Figure 4. In this situation 
heat transfer depends on whether the scale of turbulent motions is larger or smaller than 
the hot spot. Consider this situation in more detail. Turbulence is characterized by its 
injection scale L max , its dissipation scale L mi „ and its inertial range [L m ,„,L mflx ]. The 
heat transfer depends on what scales we consider the process. Figure 4 illustrates our 
point. Consider a hot spot of the size a in turbulent flow and consider Alfvenic eddies 
perpendicular to magnetic field lines. If turbulent eddies are much smaller than a, which 
is the case when a 3> L m ,„ they extend the hot spot acting in a random walk fashion. 
For eddies much larger than the hot spot, i.e. a <C L min they mostly advect hot spot. If 
a is the within the inertial range of turbulent motions, i.e. L mi „ < a < L max then a more 
complex dynamics of turbulent motions is involved. This is also the case where the 
field wandering arising from these motions is the most complex. Turbulent motions with 
the scale comparable with the hot spot induce a process of the accelerated Richardson 
diffusion (see more in §10). 

In terms of practical simulation of reconnection diffusion effects, it is important to 
keep in mind that the LV99 model predicts that the largest eddies are the most important 
for providing outflow in the reconnection zone and therefore the reconnection will not be 
substantially changed if turbulence does not have an extended inertial range. In addition, 
LV99 predicts that the effects of anomalous resistivity arising from finite numerical 
grids do not change the rate of turbulent reconnection. We note that both effects were 
successfully tested in Kowal et al. (2009). 

HEAT CONDUCTION THROUGH STREAMING OF ELECTRONS 

General considerations 

As magnetic reconnection was considered by many authors even more mysterious 
than the heat transfer in plasmas, it is not surprising that the advection of heat by turbu- 
lent eddies was not widely discussed. Instead for many year the researchers preferred to 
consider heat transfer by plasma conductivity along turbulent magnetic field lines (see 
Chandran & Cowley 1998, Malyshkin & Kulsrud 2001). This conductivity is mostly 
due to electrons streaming along magnetic field lines. Turbulent magnetic field lines al- 
low streaming electrons to diffuse perpendicular to the mean magnetic field and spread 
due to the magnetic field wandering that we discussed earlier. Therefore the descrip- 
tion of magnetic field wandering obtained in LV99 is also applicable for describing the 
processes of heat transfer. 

We start with the case of trans-Alfvenic turbulence considered by Narayan & 
Medvedev (2001, henceforth NM01). They appeal to magnetic field wandering and 
obtained estimates of thermal conductivity by electrons for the special case of turbulence 
velocity Vl at the energy injection scale L that is equal to the Alfven velocity Va- As 
we discussed earlier this special case is described by the original GS95 model and the 
Alfven Mach number Ma = (Vl/Va) = 1. We note that this case is rather restrictive, as 



the intracuster medium (ICM) is superAlfvenic, i.e. Ma > 1, while other astrophysical 
situations, e.g. solar atmosphere, are subAlfvenic, i.e. M A < 1. Different phases of 
interstellar medium (ISM) (see Draine & Lazarian 1998 and Yan, Lazarian & Draine 
2004 for lists of idealized ISM phases) present the cases of both superAlfvenic and 
subAlfvenic turbulence. 

As we discussed above, the generalization of GS95 model of turbulence for sub- 
Alfvenic case is provided in LV99. This was employed in Lazarian (2006) to describe 
heat conduction for magnetized turbulent plasmas with Ma < 1 . In addition, Lazarian 
(2006) considered heat conduction by tubulence with Ma > 1 as well as heat advection 
by turbulence and compares the efficiencies of electron heat conduction and the heat 
transfer by turbulent motions. 

Let us initially disregard the dynamics of fluid motions on diffusion, i.e. consider 
diffusion induced by particles moving along wandering turbulent magnetic field lines, 
which motions we disregard for the sake of simplicity. Magnetized turbulence with a 
dynamically important magnetic field is anisotropic with eddies elongated along (hence- 
forth denoted by ||) the direction of local magnetic field, i.e. l± < h, where _L denotes 
the direction of perpendicular to the local magnetic field. Consider isotropic injection of 
energy at the outer scale L and dissipation at the scale /j_ ;mjn . This scale corresponds to 
the minimal dimension of the turbulent eddies. 

Turbulence motions induce magnetic field divergence. It is easy to notice (LV99, 
NM01) that the separations of magnetic field lines at small scales less than the damping 
scale of turbulence, i.e. for tq < l±, m i n , are mostly influenced by the motions at the 
smallest scale. This scale l±^ m i n results in Lyapunov-type growth ~ roexp(///|| min ). This 
growth is similar to that obtained in earlier models with a single scale of turbulent 
motions (Rechester & Rosenbluth 1978, henceforth RR78, Chandran & Cowley 1998). 
Indeed, as the largest shear that causes field line divergence is due to the marginally 
damped motions at the scale around l^min the effect of larger eddies can be neglected 
and we are dealing with the case of single-scale "turbulence" described by RR78. 

The electron Larmor radius presents the minimal perpendicular scale of localization. 
Thus it is natural to associate ro with the size of the cloud of electrons of the electron 
Larmor radius r La ^ partic i e . Applying the original RR78 theory (see also Chandran & 
Cowley 1998) they found that the electrons should travel over the distance 

LRR ~ l\\,mirM{^,minl ^Lar.e) 0) 

to get separated by l^ min . 

Within the single-scale "turbulent model" which formally corresponds to Lss = 
h,min = l.L,min the distance Lrr is called Rechester- Rosenbluth distance. For the ICM 
parameters the logarithmic factor in Eq. (1) is of the order of 30, and this causes 30 
times decrease of thermal conductivity for the single-scale models 13 . 

The single-scale "turbulent model" is just a toy model to study effects of turbulent 
motions. One can use this model, however, to describe what is happening below the 



13 For the single-scale model Lrr ~ 30L and the diffusion over distance A takes Lrr/Lss steps, i.e. 
A 2 ~ LrrL, which decreases the corresponding diffusion coefficient K e , ing i e ~ A 2 / 8t by the factor of 
30. 



scale of the smallest eddies. Indeed, the shear and, correspondingly, magnetic field 
line divergence is maximal for the marginally damped eddies at the dissipation scale. 
Thus for scales less than the damping scale the action of the critically damped eddies is 
dominant. 

In view of above, the realistic multi-scale turbulence with a limited (e.g. a few 
decades) inertial range the single scale description is applicable for small scales up to 
the damping scale. The logarithmic factor stays of the same order but instead of the 
injection scale L ss for the single-scale RR model, one should use h min for the actual 
turbulence. Naturally, this addition does not affect the thermal conductivity, provided 
that the actual turbulence injection scale L is much larger than l\^ min . Indeed, for the 
electrons to diffuse isotropically they should spread from to L. Alfvenic turbulence 
operates with field lines that are sufficiently stiff, i.e. the deviation of the field lines 
from their original direction is of the order unity at scale L and less for smaller scales. 
Therefore to get separated from the initial distance of /_i_ ;nMn to a distance L (see Eq. 
(5) with Ma = 1), at which the motions get uncorrected, the electrons should diffuse the 
distance slightly larger (as field lines are not straight) than \flh. This is much larger than 
the extra travel distance ~ 30/n W( -„ originating from sub-diffusive behavior at scales less 
than the turbulence damping scale. Explicit calculations in NM01 support this intuitive 
picture. 



Diffusion for M A > 1 

Turbulence withM^ > 1 evolves along hydrodynamic isotropic Kolmogorov cascade, 
over the range of scales [L, I a] , where 

Ia^L(V a /V l ) 3 =LM a \ (2) 

is the scale at which the magnetic field gets dynamically important, i.e. V/ = Va- This 
scale plays the role of the injection scale for the GS95 turbulence, i.e. V; ~ Va(1±/ Ia) 1//3 > 
with eddies at scales less than Ia geting elongated in the direction of the local magnetic 
field. The corresponding anisotropy can be characterized by the relation between the 
semi-major axes of the eddies 

l^L(l ± /L) 2 / 3 MX l , M A >1, (3) 

where || and _L are related to the direction of the local magnetic field. In other words, for 
Ma > 1, the turbulence is still isotropic at the scales larger to Ia, but develops (1±/Ia) 1 ^ 3 
anisotropy for I < Ia- 

If particles (e.g. electrons) mean free path X 3> Ia, they stream freely over the 
distance of I a- For particles initially at distance l^min to get separated by L, the required 
travel is the random walk with the step Ia, i.e. the mean-squared displacement of a 
particle till it enters an independent large-scale eddy A 2 ~ 1^(L/Ia), where L/Ia is the 
number of steps. These steps require time 8t ~ {L/Ia)Ia/C\v c , where v parti - c / e is electron 
thermal velocity and the coefficient C\ = 1/3 accounts for ID character of motion along 



magnetic field lines. Thus the electron diffusion coefficient is 

K e = A 2 /8t^(l/3)l A v e , l A <k, (4) 

which for ^ < A constitutes a substantial reduction of diffusivity compared to its 
unmagnetized value K unmagn = Xv e . We assumed in Eq. (4) that L 3> 30/|| W( -„ (see §2.1). 

For A < < ^, K: ~ l/3 K unmagn as both the L/e/e and the additional distance 
for electron to diffuse because of magnetic field being stiff at scales less than Ia are 
negligible compared to L. For Ia — > L, when magnetic field has rigidity up to the scale 
L, it gets around 1/5 of the value in unmagnetized medium, according to NM01 . 



Diffusion for M A < 1 

It is intuitively clear that for Ma < 1 turbulence should be anisotropic from the 
injection scale L. In fact, at large scales the turbulence is expected to be weak 14 (see 
Lazarian & Vishniac 1999, henceforth LV99). Weak turbulence is characterized by 
wavepackets that do not change their , but develop structures perpendicular to magnetic 
field, i.e. decrease l± . This cannot proceed indefinitely, however. At some small scale 
the GS95 condition of critical balance, i.e. Ih/Va ~ /j_/V/, becomes satisfied. This 
perpendicular scale l trans can be obtained substituting the scaling of weak turbulence (see 
LV99) Vi ~ Vl{Ii_/L)^I 2 into the critical balance condition. This provides ltrans ~ LM^ 
and the corresponding velocity V trans ~ ViMa- For scales less than l trans the turbulence 

is strong and it follows the scalings of the GS95-type, i.e. V[ ~ Vl(L/1 1 ^)~ 1 / 3 M^/ 3 and 

/||~L(/ ± /L) 2 / 3 M~ 4/3 , M A <\. (5) 

For Ma < 1 , magnetic field wandering in the direction perpendicular to the mean 
magnetic field (along y-axis) can be described by d(y 2 )/dx ~ (y 2 )/l\\ (LV99), where 15 
Zii is expressed by Eq. (5) and one can associate lj_ with 2{y 2 ) 



3/2 

^ 2 ) 1/2 ~^ZT72 M A, h< ltrans (6) 



For weak turbulence d (y 2 ) jdx ~ LM\ (LV99) and thus 

(fy^-L^X^M 2 , I ±> ltrans- (7) 

Fig. 5 confirms the correctness of the above scaling numerically. 



The terms "weak" and "strong" turbulence are accepted in the literature, but can be confusing. As we 
discuss later at smaller scales at which the turbulent velocities decrease the turbulence becomes strong. 
The formal theory of weak turbulence is given in Galtier et al. (2000). 

15 The fact that one gets l\\ min in Eq. (1) is related to the presence of this scale in this diffusion equation. 




FIGURE 5. Root mean square separation of field lines in a simulation of inviscid MHD turbulence, 
as a function of distance parallel to the mean magnetic field, for a range of initial separations. Each 
curve represents 1600 line pairs. The simulation has been filtered to remove pseudo-Alfven modes, which 
introduce noise into the diffusion calculation. From Lazarian, Vishniac & Cho 2004. 



Eq. (6) differs by the factor M\ from that in NM01, which reflects the gradual sup- 
pression of thermal conductivity perpendicular to the mean magnetic field as the mag- 
netic field gets stronger. Physically this means that for Ma < 1 the magnetic field fluc- 
tuates around the well-defined mean direction. Therefore the diffusivity gets anisotropic 
with the diffusion coefficient parallel to the mean field K\\ part i c i e ~ \/2>K unmagn being 
larger than coefficient for diffusion perpendicular to magnetic field K± e . 

Consider the coefficient JCj_ je for Ma € 1. As NM01 showed, particles become 
uncorrelated if they are displaced over the distance L in the direction perpendicular to 
magnetic field. To do this, a particle has first to travel L RR (see Eq. (1)), where Eq. (5) 
relates /||, ra; „ and l± jTn i n . Similar to the case in §2.1, for L 3> 30h m i n , the additional 
travel arising from the logarithmic factor is negligible compared to the overall diffusion 
distance L. At larger scales electron has to diffuse ~ L in the direction parallel to 
magnetic field to cover the distance of LM\ in the direction perpendicular to magnetic 
field direction. To diffuse over a distance R with random walk of LM\ one requires 
R 2 /L 2 M^ steps. The time of the individual step is L 2 /Kj| e . Therefore the perpendicular 
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FIGURE 6. Comparison of the heat diffusion with time for hydro turbulence (left panel) and MHD 
transAlfvenic turbulence (right panel). Different curves correspond to different runs. From Cho et al. 
(2003). 



diffusion coefficient is 

K^ e = R 2 /(R 2 /[K le M 4 A ]) = K le M 4 A , M A <1, (8) 

An essential assumption there is that the particles do not trace their way back over the 
individual steps along magnetic field lines, i.e. Lrr « L. Note, that forM^ of the order 
of unity this is not accurate and one should account for the actual 3D displacement. This 
introduces the change by a factor of order unity (see above). 

TRANSFER OF HEAT THROUGH TURBULENT MOTIONS 

As we discussed above, turbulent motions themselves can induce advective transport of 
heat. Appealing to LV99 model of reconnection one can conclude that turbulence with 
M A ~ 1 should be similar to hydrodynamic turbulence, i.e. 

^dynamic ~ C dyn LV L , M A > 1, (9) 

where Cd yn ~ 0(1) is a constant, which for hydro turbulence is around 1/3 (Lesieur 
1990). This was confirmed in Cho et al. (2003) (see Figure 6 and also Cho & Lazarian 
2004) where MHD calculations were performed for transAlfvenic turbulence with M A ~ 
1. As large scale eddies of superAlfvenic turbulence are essentially hydrodynamic, the 
correspondence between the ordinary hydrodynamic heat advection and superAlfvenic 
one should only increase as M A increases. 

If we deal with heat transport, for fully ionized non-degenerate plasmas we assume 
Cdyn ~ 2/3 to account for the advective heat transport by both protons and electrons 16 . 
Thus eq. (9) covers the cases of both M A > 1 up to M A ~ 1. For M A < 1 one can 



16 This becomes clear if one uses the heat flux equation q = —K c \jT, where K c = nkBK dynamic / electr , n is 
electron number density, and kg is the Boltzmann constant, for both electron and advective heat transport. 
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FIGURE 7. (a) The textbook description of confinement of charged particles in magnetic fields; (b) 
diffusion of particles in turbulent fields; (c) advection of heat from a localized souce by eddies in MHD 
numerical simulations. From Cho & Lazarian 2004. 



estimate ^dynamic ~ d 2 C0, where d is the random walk of the field line over the wave 
period ~ to -1 . As the weak turbulence at scale L evolves over time T ~ Af^ 2 G) _1 , 
(y 2 ) is the result of the random walk with a step d, i.e. (y 2 ) ~ {x(0)d 2 . According to 
eq.(6) and (7), the field line is displaced over time % by (y 2 ) ~ LM^VaT. Combining 
the two one gets d 2 ~ LMjy L (0 1 , which provides ^ymmic ~ ^dynP^iM\, which is 
similar to the diffusivity arising from strong turbulence at scales less than l tra ns, 
K dynamic ~ c dynkransVtrans- The total diffusivity is the sum of the two, i.e. for plasma 

^dynamic ~ (P/3)LV L Mi M A < 1, (10) 

where /3 ~ 4. 



RELATIVE IMPORTANCE OF TWO PROCESSES 

General treatment 

Figure 7 illustrates the existing ideas on processes of heat conduction in astrophysical 
plasmas. They range from the heat insulation by unrealistically laminar magnetic field 
(see panel (a)), to heat diffusion in turbulent magnetic field (see panel (b)) and to heat 
advection by turbulent flows (see panel (c)). The relative efficiencies of the two latter 
processes depend on parameters of turbulent plasma. 

In thermal plasma, electrons are mostly responsible for thermal conductivity. The 
schematic of the parameter space for K part i c i e < Kdynamic is shown in Fig 8, where the 
the Mach number M s and the Alfven Mach number Ma are the variables. For Ma < 1, 
the ratio of diffusivities arising from fluid and particle motions is Kdynamicl K panicie ~ 
P<xMsMa(L/X) (see Eqs. (8) and (10)), the square root of the ratio of the electron 
to proton mass a = (m e /m p ) 1 / 2 , which provides the separation line between the two 
regions in Fig. 2, /3aM, ~ (X/L)M A . For 1 < M A < (L/X) 1 / 3 the mean free path is 
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FIGURE 8. Parameter space for particle diffusion or turbulent diffusion to dominate: application to 
heat transfer. Sonic Mach number M s is ploted against the Alfven Mach number Ma. The heat transport 
is dominated by the dynamics of turbulent eddies is above the curve (area denoted "dynamic turbulent 
transport") and by thermal conductivity of electrons is below the curve (area denoted "electron heat 
transport"). Here X is the mean free path of the electron, L is the driving scale, and a = (m e /m p ) x l 2 , 
/3 ~ 4. Example of theory application: The panel in the right upper corner of the figure illustrates heat 
transport for the parameters for a cool core Hydra cluster (point "F"), "V" corresponds to the illustrative 
model of a cluster core in Ensslin et al. (2005). Relevant parameters were used for L and X . From Lazarian 
(2006). 



less than Ia which results in K partic i e being some fraction of K, 



unmagn ■ 
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given by Eq. (9). Thus K dynamic / K partide ~ fiaM s (L/X), i.e. the ratio does not depend 

on Ma (horisontal line in Fig. 2). When Ma > (L/A) 1//3 the mean free path of electrons 
is constrained by l A . In this case K dynamic / K partide ~ fiaM s M\ (see Eqs. (9) and (4)) . 

This results in the separation line [3aM s ~ M^ 3 in Fig. 8. 



Application to ICM plasmas 



Consider plasmas in clusters of galaxies to illustrate the relative importance of two 
processes of heat transfer. Below we shall provide evidence that magnetized Intracluster 
Medium (ICM) is turbulent and therefore our considerations above should be applicable. 

It is generally believed that ICM plasma is turbulent. However, naive estimates 
of diffusivity for collisionless plasma provide numbers which may cast doubt on this 



conclusion. Indeed, in unmagnatized plasma with the ICM temperatures T ~ 10 8 K and 
and density 1(T 3 cm~ 3 the kinematic viscosity r\ unmagn ~ v ion X ion , where v ion and X ion 
are the velocity of an ion and its mean free path, respectively, would make the Reynolds 
number Re = LVL/f]unmagn of the order of 30. This is barely enough for the onset of 
turbulence. For the sake of simplicity we assume that ion mean free path coincides with 
the proton mean free path and both scale as X ~ 3T^nZ_\ kpc, where the temperature 
T3 = kT/3 keV and n_3 = n/10~ 3 cm~ 3 . This provides X of the order of 0.8-1 kpc for 
the ICM (see NM01). We shall argue that the above low estimate of Re is an artifact of 
our neglecting magnetic field. 

In general, a single value of Re uniquely characterizes hydrodynamic flows. The 
case of magnetized plasma is very different as the diffusivities of protons parallel and 
perpendicular to magnetic fields are different. The diffusion of protons perpendicular 
to the local magnetic field is usually very slow. Such a diffusion arises from proton 
scattering. Assuming the maximal scattering rate of an proton, i.e. scattering every orbit 
(the so-called Bohm diffusion limit) one gets the viscosity perpendicular to magnetic 
field ~ Vi on r Lar:ion , which is much smaller than r\ unmagn , provided that the ion Larmor 
radius r Lar ^ on <C A ;o „. For the parameters of the ICM this allows essentially inviscid fluid 
motions 17 of magnetic lines parallel to each other, e.g. Alfven motions. 

In spite of the substantial progress in understading of the ICM (see EnBlin, Vogt 
& Pfrommer 2005, henceforth EVP05, EnBlin & Vogt 2006, henceforth EV06 and 
references therein), the basic parameters of ICM turbulence are known within the factor 
of 3 at best. For instance, the estimates of injection velocity Vl varies in the literature 
from 300 km/s to 10 3 km/s, while the injection scale L varies from 20 kpc to 200 
kpc, depending whether the injection of energy by galaxy mergers or galaxy wakes is 
considered. EVP05 considers an illustrative model in which the magnetic field with the 
10 /lG fills 10% of the volume, while 90% of the volume is filled with the field of B ~ 1 
/lG. Using the latter number and assuming Vl = 10 3 km/s, L = 100 kpc, and the density 
of the hot ICM is 10~ 3 cm~ 3 , one gets Va ~ 70 km/s, i.e. Ma > 1. Using the numbers 
above, one gets Ia ~ 30 pc for the 90% of the volume of the hot ICM, which is much less 
than Xj on . The diffusivity of ICM plasma gets r\ = vu m lA which for the parameters above 
provides Re ~ 2 x 10 3 , which is enough for driving superAlfvenic turbulence at the outer 
scale L. However, as Ia increases as B 3 , Re gets around 50 for the field of 4 /iG, which 
is at the border line of exciting turbulence 18 . However, the regions with higher magnetic 
fields (e.g. 10 /lG) can support Alfvenic-type turbulence with the injection scale Ia and 
the injection velocities resulting from large-scale shear Vl(Ia/L) ~ VlM a 3 . 

For the regions of B ~ 1 /lG the value of Ia is smaller than the mean free path of 
electrons X. According to Eq. (4) the value of K e i ectr is 100 times smaller than K^ pitzer . On 



17 A regular magnetic field Bx ~ (2mkT) x l 2 c/(eX) that makes ri^/on less than X and therefore 7]x < 
Vunmagn is just 10~ 20 G. Turbulent magnetic field with many reversals over ri ar j on does not interact 
efficiently with a proton, however. As the result, the protons are not constrained until Ia gets of the order 
of r Lar,ion- This happens when the turbulent magnetic field is of the order of 2 x 10~ 9 (VL/10 3 km/s) G. At 
this point, the step for the random walk is <~ 2 x 10~ 6 pc and the Reynolds number is 5 x 10 10 . 

18 One can imagine dynamo action in which superAlfvenic turbulence generates magnetic field till l A gets 
large enough to shut down the turbulence. 



the contrary, Kj ynamic for the ICM parameters adopted will be ~ 30Ks p i tzer , which makes 
the heat transfer by turbulent motions the dominant process. This agrees well with the 
observations in Voigt & Fabian (2004). Fig. 2 shows the dominance of advective heat 
transfer for the parameters of the cool core of Hydra A ( B = 6 /iG, n = 0.056 cm 3 , 
L = 40 kpc, T = 2.7 keV according to EV06), point "F", and for the illustrative model 
in EVP05, point "V", for which 5=1^0 (see also Lazarian 2006). 

Note that our stationary model of MHD turbulence is not directly applicable to 
transient wakes behind galaxies. The ratio of the damping times of the hydro turbulence 
and the time of straightening of the magnetic field lines is ~ . Thus, for Ma > 1, 
the magnetic field at scales larger than Ia will be straightening gradually after the hydro 
turbulence has faded away over time L/Vl. The process can be characterized as injection 
of turbulence at velocity Va but at scales that increase linearly with time, i.e. as Ia + Va?. 
The study of heat transfer in transient turbulence and magnetic field "regularly" stretched 
by passing galaxies is an interesting process that requires further investigation. 



RICHARDSON DIFFUSION AND SUPERDIFFUSION ON SMALL 

SCALES 



All the discussion above assumed that we deal with diffusion within magnetized 
plasmas over the scales much larger than the turbulence injection scale L. Below we 
show that on the scales less than L we deal with non- stationary processes. 



The advection of heat on scales less than the turbulent injection scale L happens 
through smaller scale eddies. Thus the earlier estimate of turbulent diffusion of heat in 
terms of the injection velocity and the injection scale does not apply. In the lab system of 
reference the transfer of heat is difficult to describe and one should use the Lagrangian 
description. 

One can consider two-particle turbulent diffusion or Richardson diffusion by dealing 
with the separation £{t) = x(t) — x'(t) between a pair of Lagrangian fluid particles (see 
Eyink et al. 2011). It was proposed by Richardson (1926) that this separation grows in 
turbulent flow according to the formula 



with a scale-dependent eddy-diffusivity Kdynamici^) ■ I n hydrodynamic turbulence 
Richardson deduced that JQfy Maw( - c (£) ~ £ 1//3 ^ 4 / 3 (see Obukhov 1941) and thus 
£ 2 (t) ~ £t 3 . An analytical formula for the 2-particle eddy-diffusivity was derived 
by Batchelor (1950) and Kraichnan (1966): 



Richardson-type advection of heat 




(ii) 




(12) 



with 8Ui(£,t) = Ui(x + £,t) — Ui(x,t) the relative velocity at time t of a pair of fluid 
particles which were at positions x and x + 1 at time 0. 

How can one understand these results? Consider a hot spot of the size / in a turbulent 
flow. The spot is going to be mostly expanded by turbulent eddies of size /. The turbulent 
velocity u(l) = jjjl(t) for Kolmogorov turbulence is proportional to / 1//3 . Performing 
formal integration one gets an asymptotic solution for large time scales l 2 (t) ~ ? 3 , which 
corresponds to the Richardson diffusion law. Physically, as the hot spot extends, it is 
getting sheared by larger and eddies, which induce the accelerated expansion of the hot 
spot. 

For magnetic turbulence the Kolmogorov-like description is valid for motions in- 
duced by strong Alfvenic turbulence in the direction perpendicular to the direction of 
the local magnetic field 19 . Thus we expect that Richardson diffusion to be applicable to 
the magnetized turbulence case. 



Superdiffusion of heat perpendicular to mean magnetic field 

The effects related to the diffusion of heat via electron streaming along magnetic field 
lines are different when the problem is considered at scales 3> L and <C L. This difference 
is easy to understand as on small scales magnetized eddies are very elongated, which 
means that the magnetic field lines are nearly parallel. However, as electrons diffuse into 
larger eddies, the dispersion of the magnetic field lines in these eddies gets bigger and 
the diffusion perpendicular to the mean magnetic field increases 20 

SuperAlfvenic turbulence: 
On scales k^ 1 < Ia, i.e., on scales at which magnetic fields are strong enough to influence 

turbulent motions, the mean deviation of a field in a distance k^ 1 = 8z is given by LV99 
as 

< W > = 3 3/2 L l/2 ' Ma > 1 < 13) 

Thus, for scales much less than L (see also Yan & Lazarian 2008) 

Ke,±™(j£j k^^^^Ak^^k^Ia)- 1 , M a >1, (14) 

which illustrates the non-stationary regime of superdiffusion, where the diffusion coef- 
ficient changes with the scale k~l . 

Sub Alfvenic turbulence: 
On scales larger than l tr , the turbulence is weak. The mean deviation of a field in a 



ly The local magnetic field direction fluctuates in the lab system of reference. Thus the results of the 
diffusion in the lab system are less anisotropic. 

20 Below we consider turbulent scales that are larger than the electron mean free path X e . Heat transfer 
at smaller scale is not a diffusive process, but happens at the maximal rate determined by the particle 
flux nv t k provided that we deal with scales smaller than The perpendicular to magnetic field flux is 
determined by the field line deviations on the given scale as we discussed above (see also LV99). 



distance Sz is given by Lazarian (2006): 



< (8x) 2 > 1/2 = ^^Ml M A <\. (15) 



For the scales L > k,, 1 = 8z we combine Eq. (15) with 



8z = ^kappa e ^8t (16) 
and get for scales much less than L 



8x 2 K e nSz , , , , a 

Ke,± ~ -fr = ~^r M A ~ ^n^iL)-^, (17) 



which for a limiting case of k e \\ ~ L coincides up to a factor with the Eq. (8). 

Eqs. (14) and (17) certify that the perpendicular diffusion at scales much less than 
the injection scale accelerates as z grows. 



Comparison of processes 

Both processes of heat transport at the scales less than the turbulence injection scale 
are different from the diffusion at large scales as the rate of transport depends on the 
scale. However, the description of heat transport by electrons is more related to the mea- 
surements in the lab system. This follows from the fact that the dynamics of magnetic 
field lines is not important for the process and it is electrons which stream along wan- 
dering magnetic field lines. Each of these wandering magnetic field lines are snapshot 
of the magnetic field line dynamics as it changes through magnetic reconnection its 
connectivity in the ambient plasma. Therefore the description of heat transfer is well 
connected to the lab system of reference. On the contrary, the advection of heat through 
the Richardson diffusion is a process that is related to the Lagrangian description of the 
fluid. Due to this difference the direct comparison of the efficiency of processes is not 
so straightforward. 

For example, if one introduces a localized hot spot, electron transport would produce 
heating of the adjacent material along the expanding cone of magnetic field lines, while 
the turbulent advection would not only spread the hot spot, but also advect it by the 
action of the largest eddies. 



OUTLOOK ON THE CONSEQUENCES 

Magnetic thermal insulation is a very popular concept in astrophysical literature dealing 
with magnetized plasmas. Our discussion above shows that in many cases this insulation 
is very leaky. This happens due to ubiquitous astrophysical turbulence which induces 
magnetic field wandering and interchange of pieces of magnetized plasma enabled by 



turbulent motions. Both processes are very closely related to the process of fast magnetic 
reconnection of turbulent magnetic field (LV99). 

As a result, instead of an impenetrable wall of laminar ordered magnetic field lines, 
the actual turbulent field lines present a complex network of tunnels along which elec- 
trons can carry heat. As a result, the decrease of heat conduction amounts to a factor in 
the range of 1/3 for mildly superAlfvenic turbulence to a factor ~ 1/5 for transAlfvenic 
turbulence. The cases when heat conductivity by electrons may be suppressed to much 
greater degree include highly superAlfvenic turbulence and highly subAlfvenic turbu- 
lence. In addition, turbulent motions induce heat advection which is similar to turbulent 
diffusivity of unmagnetized fluids. 

The importance of magnetic reconnection cannot be stressed enough in relation to 
the process of heat transfer in magnetized plasmas. As a consequence of fast magnetic 
reconnection plasma does not stay entrained on the same magnetic field lines, as it is 
usually presented in textbooks. On the contrary, magnetic field lines constantly change 
their connectivity and plasma constantly samples newly formed magnetic field lines 
enabling efficient diffusion. Therefore we claim that the advection of heat by turbulence 
is an example of a more general process of reconnection diffusion. It can be noticed 
parenthetically that the turbulent advection of heat is a well knows process. However, for 
decades the discussion of the process avoided in astrophysical literature due the worries 
of the effect of reconnection that inevitably should accompany it. The situation has 
changed with better understanding of magnetic reconnection in turbulent environments 
(LV99). It worth pointing out that our estimates indicate that in many astrophysicaly 
important cases, e.g. for ICM, the advective heat transport by dynamic turbulent eddies 
dominates thermal conductivity. 

Having the above processes in hand, one can describe heat transport within magne- 
tized astrophysical plasmas. For instance, we discussed the heat transfer by particle and 
turbulent motions for Ma < 1 and Ma > 1 . It is important that we find that turbulence 
can both enhance diffusion and suppress it. We showed that when X gets larger than Ia 
the conductivity of the medium ~ M A 3 and therefore the turbulence inhibits heat trans- 
fer, provided that K e > K^ynamic- Along with the plasma effects that we mention below, 
this effect can, indeed, support sharp temperature gradients in hot plasmas with weak 
magnetic field. 

As discussed above, rarefied plasma, e.g. ICM plasma, has large viscosity for mo- 
tions parallel to magnetic field and marginal viscosity for motions that induce perpen- 
dicular mixing. Thus fast dissipation of sound waves in the ICM does not contradict the 
medium being turbulent. The later may be important for the heating of central regions of 
clusters caused by the AGN feedback (see Churasov et al. 2001, Nusser, Silk & Babul 
2006 and more references in EV06). Note, that models that include both heat transfer 
from the outer hot regions and an additional heating from the AGN feedback look rather 
promissing (see Ruszkowkski & Begelman 2002, Piffaretti & Kaastra 2006). We predict 
that the viscosity for 1 /iG regions is less than for 10 /lG regions and therefore heating 
by sound waves (see Fabian et al. 2005) could be more efficient for the latter. Note, that 
the plasma instabilities in collisionless magnetized ICM arising from compressive mo- 
tions (see Schekochihin & Cowley 2006, Lazarian & Beresnyak 2006) can resonantly 
scatter particles and decrease X . This decreases further K e compared to K unmagn but in- 



creases Re. In addition, we disregarded mirror effects that can reflect electrons back 
(see Malyshkin & Kulsrud 2001 and references therein), which can further decrease K e . 
While there are many instabilities that are described in plasmas with temperature gradi- 
ent, many of those are of academic interest, as they do not take into account the existence 
of ambient turbulence. 

For years the attempts to describe heat transfer in magnetized plasma were focused 
on finding the magic number which would be the reduction factor characterizing the 
effect of magnetic field on plasmas' diffusivity. Our study reveals a different and more 
complex picture. The heat transfer depends on sonic and Alfven Mach numbers of 
turbulence and the corresponding diffusion coefficient vary substantially for plasmas 
with different level of magnetization and turbulent excitation. In different astrophysical 
environments turbulence can both inhibit or enhance diffusivity depending on the plasma 
magnetization and turbulence driving. 

The issues of "subdiffusivity" or magnetic field retracing their paths was a worrisome 
issue that for years impeded the progress in understanding heat transport in plasmas. We 
claim that the retracing does happen, but on the scales which are of the order of the 
eddies at the dissipation scale. As an electron has a finite Larmor radius in the retracing 
the same magnetic field line it experiences the deviations from its original trajectory. On 
the scale less than the dissipation scale these deviations grow from the electron Larmor 
radius in accordance with Lyapunov exponents, but on larger scale the separation is 
determined by field wandering only and does not depend on the Larmor radius. Thus 
the effect of retracing for heat transfer in real-world astrophysical turbulence with a 
substantial separation of the turbulence injection scale and dissipation scales is marginal. 

On the contrary, the issue of "superdiffusivity" may be important for heat transfer 
on the scales less than the turbulence injection scale. Richardson diffusion or more 
correctly its anisotropic analog present in magnetized plasma (see Eyink et al. 2011) 
is an example of superdiffusion induced by eddies of increasing size. A similar effect is 
also true for magnetic field line wandering. The effect of "superdiffusive" heat transfer 
requires additional studies. 

It is worth mentioning that another parameter that determines the heat flux into the 
magnetized volume is the area of the contact of plasmas with different temperatures. 
For instance, if the magnetic flux is "shredded", i.e. consists of numerous separated 
individual flux tubes, then the heating of plasma within magnetized tubes may be more 
efficient. For instance, Fabian et al. (201 1) appealed to reconnection diffusion of ambient 
plasma into "shredded" magnetic flux of NGC1275 in Perseus cluster in order to explain 
heating and ionization of the magnetic filaments. 

In view of the discussion above one can conclude that realistically turbulent mag- 
netic fields do not completely suppress heat conductivity of astrophysical plasmas. The 
decrease of thermal conductivity depends on the Alfven Mach number of turbulence. At 
the same time, turbulent motions enhance heat transport via heat advection. In special 
situations, e.g. in very weakly turbulent magnetic field, the transport of heat in plasmas 
may still be slow. 



21 Many of these papers do not use realistic models of turbulence and therefore overestimate the effect of 
electron reflection. 
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